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Abstract: Recent data from cosmic ray experiments such as PAMELA, Fermi, ATIC and 
PPB-BETS all suggest the need for a new primary source of electrons and positrons at high 
(> 100 GeV) energies. Many proposals have been put forth to explain these data, usually 
relying on a single particle to annihilate or decay to produce . In this paper, we consider 
models with multiple species of WIMPs with significantly different masses. We show if such 
dark matter candidates Xi annihilate into light bosons, they naturally produce equal anni- 
hilation rates, even as the available numbers of pairs for annihilation n^. differ by orders of 
magnitude. We argue that a consequence of these models can be to add additional signal 
naturally at lower (~ 100 GeV) versus higher (~ TeV) energies, changing the expected spec- 
trum and even adding bumps at lower energies, which may alleviate some of the tension in 
the required annihilation rates between PAMELA and Fermi. These spectral changes may 
yield observable consequences in the microwave Haze signal observed at the upcoming Planck 
satellite. Such a model can connect to other observable signals such as DAMA and INTE- 
GRAL by having the lighter (heavier) state be a pseudo-Dirac fermion with splitting 100 keV 
(1 MeV). We show that variations in the halo velocity dispersion can alleviate constraints 
from final state radiation in the galactic center and galactic ridge. If the lighter WIMP has a 
large self-interaction cross section, the light- WIMP halo might collapse, dramatically altering 
expectations for direct and indirect detection signatures. 
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1. Introduction 

Signals from a wide range of sources have prompted a rethinking of the nature of dark matter. 
The positron fraction of cosmic rays has shown a steep rise |l]], implying the presence of a 
new primary source of antimatter. The Fermi ^ 3 ATIC |J and PPB-BETS j| experiments 
all show a harder spectrum of e^at high 100 GeV) energies than had been predicted from 
secondary production, with a break at ~ 1 TeV seen at Fermi and HESS || ||. Data from 
INTEGRAL suggest a new source of low energy positrons, while the DAMA modulation 
signal has extended to 8.2a. 

Remarkably, there is a simple "unified" framework that explains all of these anomalies (?j . 
One of the principle changes is that rather than annihilating to standard model channels, the 
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WIMP can annihilate dominantly into a new light, unstable boson, which in turn maintains 
equilibrium with the standard model and decays through a small mixing parameter ||||. The 
light mass kinematically prevents production of antiprotons § , allowing consistency with the 
PAMELA antiproton measurements [fLG[ |. Such scenarios have been widely discussed in the 
context of present dark matter anomalies || |ll|, ||, |7], [l2|, 13, 14, [l5| . 



In SUSY models, the LSP is stable by R-parity, and thus there is a unique WIMP 
candidate. In contrast, in these models, the symmetry stabilizing the dark matter against 
decay is not R-parity, and thus there is no compelling reason that there should be a single 
stable neutral particle. Indeed, it is very possible - perhaps even likely - that there are a 
variety of parities in the theory, and that many different particles together comprise the dark 
matter. We shall see that in models where dark matter freezes out by annihilation into a 
light boson, the annihilation rates of different components of dark matter are naturally equal, 
even though r? x can differ between the species by orders of magnitude. As a consequence, 
the spectra of positrons and electrons can have interesting structures, including crests and 
troughs. This phenomenon arises naturally, without any tuning of the parameters. Such a 
model can alleviate some tension between PAMELA, which seems to require somewhat larger 
cross sections at lower energy, and Fermi, which requires lower cross sections but extending 
to higher energy. 

In such models, one can explain the INTEGRAL signal with a heavy (~ TeV) mass 
state exciting into an MeV excited state as in ||, while explaining DAMA through a lighter 
(~ 100 GeV) state exciting into a ~ 100 keV excited state, as in |E6|, [L7|]. This scenario, to 
employ multiple states for inelastic and exciting dark matter (which we refer to as MiXDM) 



was proposed in IS]. Although DAMA and INTEGRAL both play roles in motivating the 
spectrum considered, we shall not focus on these aspects, focusing instead on the higher 
energy cosmic ray signals (and their synchrotron counterparts). Because of the motivation in 
[jnj, we shall generically refer to these models, in which multiple states annihilate comparably 
into light bosons, as "MiXDM" models. 

In this paper, we shall explore the possibility of testing such a scenario. In section ^, 
we will review the setup and in section || the model building issues, including models where 
different annihilation channels arise. In section [2], we discuss the cosmic ray signatures from 
such a scenario, and discuss the interesting spectral features that can arise. In [5], we go 
beyond the e + e~ signals, we shall argue that signatures in the microwave haze of the inner 
galaxy may help to show if the dark sector has multiple components. In section |6| we consider 
further generalizations, with more states, while in section 0, we note the important changes 
to direct detection experiments that can arise. Finally, we conclude in section ||. 

2. Dark Matter with Multiple Components 



The idea that dark matter has multiple components is not new. In [19, 2C], it was argued 
that a subdominant component of dark matter in the halo could still give significant signals 
at direct detection experiments, and indirect signals related to WIMP capture (such as in the 
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Earth). fl2l[ argued that a sector combining a TeV WIMP together with a GeV mass WIMP 



could explain PAMELA and DAMA. As we have already noted, [18] constructed a model in 
which a TeV particle was responsible for the high energy signatures, as well as the 511 keV 
emission observed by INTEGRAL, while a separate particle with a 100 keV splitting would 
explain DAMA. 

Our effort here is not to focus on the direct detection signals of multiple states, but rather 
their indirect signals. We shall see that if dark matter freezes out and annihilates presently 
into a light force carrier, the annihilation rates are naturally equal. Because the lighter WIMP 
has its signal compressed into a narrower energy range, this can lead to an enhancement of 
e^at lower energies, possibly even yielding a bump in the e + or e + + e~ spectrum. 

2.1 Annihilation rates 





Figure 1: Perturbative annihilation of dark matter into <f>. 

Let us begin by sketching out the process described in [j^] whereby dark matter freezes 
out through annihilations into a new metastable force carrier cp, to which it couples with 
strength g. The annihilation proceeds through the diagram shown in figure [l], and has a cross 
section (at high energy) that scales as 



av ~ 



a 

Mj 



(2.1) 



The dark matter stays in thermal equilibrium with the standard model through the interac- 
tions of 4> with the standard model, in particular the process e<j) «-> e^f or (jxp <-> ff, depending 
on the model. As was emphasized in m, O, 22, E3fl, this allows the "WIMP miracle" to func- 



tion, even though the WIMP has no strong interactions to the standard model. Thus, one 
has the usual relation [24] that the dark matter density scales inversely with cross section. 
I.e., 



(av) 



Ml 



(2.2) 
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where we focus only on the dominant mass scaling. The annihilation rate of the dark matter 
in the present universe then simply scales as 



n * ^ ~ ( 1!% ~ constant - ( 2 - 3 ) 

Remarkably, because the number density scales as mass, the leading dependencies cancel 
precisely. Thus one sees immediately that in scenarios where the dark matter freezes out 
through annihilations into a light force carrier, the annihilation rates of species of different 
mass are naturally equal. 1 

There are of course caveats to this. We are assuming that the coupling g is the same 
between the states. This arises automatically if is a boson of a new gauge group under 
which Xi have the same representation, or if there is a symmetry relating the different Xii 



broken softly by their masses. The Sommerfeld effect |25|, [26|, 27, 28, 23, 12 1, which is often 
invoked to explain the high cross sections needed to explain the cosmic ray data, can for 
instance have already saturated for one WIMP and not the other, or be in a v~ 2 resonance 
region for one and not the other. Factors of two can arise as well for fermions versus scalars, 
or if the representations of the groups are not identical. Finally, there is the usual logarithmic 
dependence on mass. That said, it is important to note that the starting point for considera- 
tion is that the annihilation rates should be comparable, and such effects can be invoked to 
move away from that. 



3. Models of MiXDM 



Scenarios where multiple states are present are trivial to construct. Probably the easiest 
means to do this is to WIMPs which are charged under a new U(l)darki with gauge field </>„, 
as was first considered by p9|| . To this we add a small < GeV mass to the gauge boson, and 
a small splitting 5 between the Majorana states for the WIMPs. 

£ = Xi PXi + \^ u F d ^ + eF^F*" + m 2 ^f + MmXi + **XiXi- (3-1) 
One can also trivially replace the fermionic WIMP with a complex scalar. 

C = {D^iTD^Xi + \Fi v F d ^ + eF^F*" + m 2 ^ + M 2 iX *Xi + M^XiXi- (3.2) 

Here i indexes the WIMP species, and we have suppressed the necessary dark Higgs mecha- 
nism. Reduced to a single species with 5 ~ lMeV, this is just the original XDM model pro- 
posed by ||. To explain INTEGRAL, such a model requires a double excitation. If the Milky 
Way halo has a uniform and low velocity dispersion, it is difficult to achieve this, as pointed 

1 Note that while we are focused on models with new force carriers, this would also be true for annihilations 
into W's, Z's, or even light fermions, so long as the WIMP mass were suitably large compared to the particle 
into which it annihilated. 
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out by |30|], prompting development of single excitation models ]|, |?]] and inverted models 
[fjl| , |32]| , which more easily satisfy the INTEGRAL rate due to the naturally large scattering 
cross section from the light mediator M. Some recent simulations involving baryons show a 



significant increase of the velocity dispersion in the galactic center [33, [54|, 35, ^g, 37], which 
can alleviate these constraints |38| , |39|| , possibly opening the possibility of double-excitations, 
although this is unclear. 

Meanwhile, the interaction with the standard model can allow us to construct models 
of inelastic dark matter Since heavier WIMPs have lower modulation fraction, the con- 



straints on them are stronger [17, 40, 41, 42, 43 1. While simulations and simulation- inspired 



halo models allow for heavier WIMPs 4C , 44 ] , having two WIMPs is a simple approach to this 
problem as well. Note that the annihilation cross section goes as o:^ ark /M? while the scat- 
tering cross section scales as e 2 /m1. Thus, even a subdominant WIMP can have significant 



scattering cross sections, even beyond the natural "order one" level considered by fl9|| . 

Thus, taking the above Lagrangian with two species one arrives at the MiXDM scenario, 
which, suitably supersymmetrized, is just the setup proposed by |i~q] . 

We note that if there is an approximate symmetry protecting the WIMPs from both M 
and 5, it is natural that 5«/Mj is constant, or at least roughly so. 

An interesting generalization of this would be to a situation with two U(l)'s in the dark 
sector (or even more generally, each WIMP would have its own U(l)). I.e., we would modify 
by 



(3.3) 



In this case, it is reasonable to expect that the soft breaking of the masses for x would 
radiatively feed into determining the masses of the 4>. Thus, one would expect rrii/Mi = 
constant. This could be generated for instance when the WIMPs, themselves, serve as gauge 
mediators for the masses of the 



's as described by p3| , 45], but not in cases where the 
mediation was dominated by the effective D-terms, as in [46, |l^, 47]. 

Finally, we can consider situations where ^ is a scalar, and there is a Z2 symmetry 
between the two WIMPs. I.e., we could consider 



£ = Xi 



d^d^cf) + \4> 2 h)h + m 2 (t> 2 



+ Mi)xx + $iXX 



(3.4) 



Which can be generalized to a multi-(/> theory trivially, in which case again, one could argue 
that the <pi masses would be proportional to the Xi masses. The scalar mediated theories 
can yield XDM phenomenology, but is difficult to produce the iDM phenomenology naturally 
pg| ] . Since our focus is on cosmic rays, the particular Lorentz representation is not especially 
crucial, but a scenario with the DAMA signal likely fits best with vectors. 

We should note that there are a wide variety of studies on the possibility of light force 
carriers. [49] first studied them in the context of MeV scale dark matter, as an explanation 
of INTEGRAL. While XDM is limited to ~ TeV masses for X and 100 MeV ~ GeV for c/>, 
[|i~l| ] considered the setup for a much wider kinematical range, including MeV WIMPs and 
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~ 20 GeV force carriers. Decaying models have also invoked new light force carriers |5l|] 
While we focus on cosmic ray signals here, a broad set of collider and other accelerator studies 

1 will be 



as well as astrophysical signatures |58|, 59, 60, 61, 62 



essential in testing this picture. (For a broad review of existing and future limits, see |60| .) 



4. Electron and Positron Signals of Multiple WIMPs 

Before progressing to consider in details the electron signals from multiple WIMPs, we can 
sketch out what sort of range of signals is possible. The simplest scenario is to consider two 
WIMPs, of mass M\ ~ 100 GeV and M<i ~ 1 TeV, annihilating through the channel XX ~^ ^ 
with <j) — > e + e~ (which naturally occurs for < 2m^). We calculate the cosmic ray spectra 
as described in appendix [A|, and show the results in figure |2[ Here we see immediately 
the important consequence of this scenario. Although the standard case of a heavy WIMP 
annihilating through this channel yields very hard spectra, the inclusion of the lighter WIMP 
leads to a softening. Indeed, with a slightly higher rate than equal (which can arise, as 
described, from scalar /fermion differences, or differences in the Sommerfeld effect), bump like 
features can appear in the spectrum. This is possibly the most important point of this section: 
the presence of bumps or a flattening in the positron spectrum may indicate the presence of a 
second, lighter WIMP. Although no such features have been found at this point, future data 
releases may show them as PAMELA goes to higher energies and with the higher statistics 
and energies of AMS-02 @. 2 

(a) (b) 




Energy (GeV) Energy (GeV) 

Figure 2: Fermi/HESS electron plus positron (left) and PAMELA positron (right) signals from 
annihilation of 2 XDM species with M Xl = 120 GeV and M X2 = 1.2 TeV through e ± channel. 
BF specifies the cross section of the heavier WIMP while x : y specifies the lighter WIMP has an 
annihilation rate x/y times that of the heavier. 

2 Note that the parameters chosen here, as in the examples through this paper, are consistent with the CMB 
limits described in E^, [57J . 
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Such "bumpy" features have been considered before in the context of pulsars |68|, 69. I? 



where bumps are associated with the duration between the initial release of positrons from 
the pulsar wind nebula and the present time. Although it is likely difficult to get some of 
the truly "spikey" signals considered by |6^, 39, [70], it is intriguing that broader features can 
be realized quite naturally for annihilating WIMP models, and are even expected in theories 
with multiple WIMPs (see additionally |7l|] for a discussion of such features from multi-WIMP 
setups in decaying models) as well as pulsars. 

While in figure || we have varied the relative annihilation rates slightly, we can also vary 
the mass of the lighter state, which we do in figure ^. In cases where the lighter WIMP has 
a slightly larger cross section (plotted 2:1) the "feature" produced in the e ± spectrum would 
seem in tension with the Fermi data if it is too high (^ 150 GeV). Lighter WIMPs seem to 
have a good fit with the Fermi data, and it is intriguing to speculate on the change in slope 
in the experiment Fermi electron data. This slope change, while not statistically significant, 
could arise from the break associated with a ~ 100 GeV WIMP into light force carriers. 
Should this persist in the future, it may be pointing to the MiXDM scenario. 



(a) 



(b) 



XDM eV Channel, with M Zi =1 00 or 1 20 or 1 50 or 200 GeV 
and IVL =1 .2 TeV at 2:1 
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Figure 3: Cosmic ray signals as in Fig. | with M X1 
through the XDM e ± channel. 



100, 120, 150, 200 GeV and M x 



1.2 TeV 



4.1 Annihilations with Different Force Carriers 

As noted above, it is possible to have different WIMPs, each with its own </>j. It is natural to 
expect in this case that the <fii would receive different radiative corrections to their masses, 
which could naturally yield a situation where Mj : <5j : rrii was constant. 

As a means to study this, we will consider the situation where the lighter particle anni- 
hilates as 2xi — ► 20i and then the mediator decays to leptons (pi — > , while the heavier 
annihilates to 2\2 — * 2<p2 where <p2 24>2 that 4>' 2 then decays to . For the analytical 
formulas for the injection spectra of the one can seej72], ff^ j. 
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In Fig. |] we show the results where I = e or I = /i. This results in a very smooth 
change of the e + + e~ flux at E = M X1 relative to that of just xii f° r cases where relative 
annihilation rates of 1:1 are assumed. Both cases provide good fit to the Fermi data while 
their projected (B.F. calculated from fitting only to the Fermi data) positron fraction agrees 
with the PAMELA data. 

In Fig. (|a,b we show the case where M Xl = 200GeV annihilating to 2(f>i with < Im^, 
thus 4>\ — ► e + e~ (XDM e + e~ channel) and M X2 = 1.6TeV annihilating to 202 where 02 — * 20 2 
with m^ 2 ~ Wm^ and (f)^ — ► e + e~ as m^> < 2to^(XDM 2-step e + e~ channel). In Fig. ||c,d 
we study the case where M X1 = 300GeV and M X2 = 3.6TeV with Im^ < m^^m^ < 2m- K , 
thus the final product of the annihilations of the two species are (4>i,4>' 2 — ► fA + fj,~), 

(XDM channel for the lighter and XDM 2-step channel for thr heavier species). 

Finally in Fig. ||e,f we study the case where M Xl = 250GeV annihilating through XDM e + e~ 
channel and M X2 = 2.7TeV through XDM fi + channel. Our choices for M X2 are suggested 
by the need to fit to the Fermi data and not violate the HESS upper limits on the e + e~ flux. 
As 2-body cascade decays of (ft's and 3-body decays of muons produce significantly softer 
injection spectra of than XDM e + e~ , higher masses of M X2 are needed to fit to the Fermi 
data. 

In viewing the various examples of MiXDM annihilation signals, there are a few important 
issues that present themselves: first, while 4:1 rates typically clearly exhibit a feature in the 
electron and/or positron spectra, they are also generally in strong tension with the current 
Fermi data as well. Signals with 2:1 and smaller ratios are more difficult to distinguish in 
the electronic spectra. At the same time, the projected positron fraction fits the PAMELA 
data well for ratios of 1:1 and 2:1, alleviating some existing tension between the suggested 
annihilation rates from Fermi vs PAMELA data in single WIMP scenarios. Finally, should 
such a signal be present in the upcoming experimental results, it will still be unclear whether 
this is arising from dark matter, or from a conventional astrophysical source, such as pulsars. 
To answer the question whether a second WIMP is present in the data, we seek an additional 
means to constrain the signal. 



4.2 Inverse Compton Signals at Fermi 



As has been discussed elsewhere [74, 75, 76, 77, 78, 79, 80], models of WIMPs that explain 



PAMEL A/Fermi also yield a substantial gamma ray signal in the galactic center region from 
inverse-Compton scattering. Remarkably, based upon the Fermi data release [^lj have argued 
that such a signal is present. Because MiXDM models have an additional soft component, 
this can modify the spectrum of gamma rays. 

We show an example of the gamma ray spectrum for a MiXDM model in Fig ||. A 
somewhat smaller boost in the Haze region seems to give an additional contribution to the 



gamma rays which is at high energies comparable to what is extracted in [81]. Such a variation 
could be easily accounted for by a v" 1 scaling, as shown. 
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and XDM 2-step to eV Channel, with WL =1 .6 TeV 








"'"Sou. 1 -,- 


HESS syst. uncert. in energy 




HESS Data 


"""--„. '\ 


i — i — i Fermi Data 




BF Xj = 120, 4:1 




BF^"- 140, 2:1 .p> 




BF^ = 160, 1:1 J^s'' 




BF X " = 160, V2s ' Z-^p??-''' 




Background/'/.- \-'"' 




/ / // DM Contribution \ \\ 

///.■•' 1 1 '» 



BF Xi = 120, 4:1 
BF =140, 2:1 

. BF X ~ = 160, 1:1 
BF = 160, 1:2 
Background 

1 PAMELA Data 




XDM eV Channel, 
with M % =200 GeV 

and XDM 2-step e + e" Channel, 
with M„ =1.6 TeV 



Energy (GeV) 
(c) 

XDM U.V Channel, with M z =300 GeV 

and XDM 2-step to nV Channel, with M Xi =3.6 TeV 



HESS syst. uncert. 
i HESS Data 
i Fermi Data 

= 920, 4:1 
= 1100, 2:1 
= 1100, 1:1 
= 1200, 1:2 
Backgrounds-^ 





XDM eV Channel, with M % =250 GeV 
and XDM Channel, with M =2.7 TeV 




HESS syst. uncert 
i HESS Data 
h Fermi Data 

= 310, 4:1 
s 470, 2:1 
= 610, 1:1 



Background 



Energy (GeV) 
(d) 



: 920, 4:1 
.1100, 2:1 
: 1100, 1:1 
: 1200, 1:2 
Background 
1 PAMELA Data 




XDM nV Channel, 
with M =300 GeV 

and XDM 2-step u. + (.i Channel, 
with M.. =3.6 TeV 



Energy (GeV) 
(f) 



BF Tn = 310, 4:1 
^ = 470, 2:1 
= 610, 1:1 
= 690, 1 :2 

Background 

i— •— i PAMELA Data 




XDM eV Channel, 
with M =250 GeV 

and XDM u\ + n~ Channel, 
with M v =2.7 TeV 



Energy (GeV) 



Energy (GeV) 



Figure 4: Cosmic ray signals as in Fig. g. Top: XDM e ± with M Xl = 200 GeV and XDM 2-step 
to e ± with M X2 = 1.6 TeV. Middle: XDM /i ± with M Xl = 300 GeV and XDM 2-step to ^ ± with 
m X2 = 3.6 TeV. Bottom: XDM e ± with m Xl = 250 GeV and XDM // ± with m X2 = 2.7 TeV. 



4.3 General fits and constraints 

In addition to the examples presented, we have also studied a variety of combinations of XDM 
e ± , XDM /i ± , XDM n ± 7 XDM 2-step to e ± , XDM 2-step to ^ and combinations of XDM to 
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Figure 5: A comparison of the residual Fermi Haze data \l\ < 10° and 10° < \b\ < 30° from |pl[ , 
versus DM ICS contributions for uniform boost and an Einasto profile (left) and with a u _1 scaling 
assuming Vdisp oc r -1 / 4 (right). 



and with various annihilation ratios. In general, for each of these channels we could 
find a combination of masses Mi and M2, that fit the Fermi e + + e~ flux data within 95% 
C.L. Thus - and not surprisingly - such combinations of two annihilating species fit the Fermi 
e + + e~ data to 95% C.L. in a fairly generic manner. 

While we have demonstrated good agreement between the Fermi/PAMELA data sets 
from MiXDM models, we should consider other constraints as well. A variety of constraints 
from the diffuse gamma-ray background p^ , |83| ] can arise, but the broad uncertainties in halo 
models allow the models discussed here to be unconstrained. In contrast, there are limits 
from the linear regime of the universe and the local halo that more directly connect to the 
local observations, or do so with fewer uncertainties about galaxy formation. 



The most model-independent of these come from the CMB. Recently, |65j, 66, S7 have 
shown the CMB can constrain the maximum annihilation into e + e~ in the early universe. At 
z ~ 10 3 those electrons will heat and ionize the photon-baryon plasma, changing the power 
spectrum and the polarization of the CMB. Thus constraints can be put on the annihilation 
rate at z ~ 10 3 based on the WMAP5 data. 

Considering that the velocity dispersion of the DM at z ~ 10 3 is lower than what we 
expect locally, those constraints on the DM annihilation can be considered as the saturated 



values of the Sommerfeld enhancement (although note |84[]). [37j have calculated the maxi- 
mum values of allowed DM annihilations for the XDM e^, XDM fj^ and XDM tt^ channels. 
In all the cases presented in Fig. ^ to [| the suggested B.F. are below or similar to the max- 
imum allowed Sommerfeld enhancement allowed by the WMAP5 data. Although there are 
some additional uncertainties in relating the local signal to that of the early universe (in 
particular, the local density and contributions from any local substructure), future data from 
Planck should constrain further the various annihilation channels [67]. 
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Figure 6: A comparison of the signal from FSR+ICS versus HESS data in the Galactic Ridge, 
normalized to local PAMELA and Fermi signals assuming a cross section scaling as v°, v" 1 and v~ 2 
with Vdisp oc r^ 1 / 4 , using an Einasto profile with a = 0.17. The lowest data point from Fermi is 
extracted from the maps of [pl|. 



Dark matter annihilation into charged particles naturally produces high energy 7's from 
final state radiation (FSR) which can place strong constraints in the galactic center (GC) and 
galactic ridge (GR) (8^, |86|, |87j. These limits have been shown to be generally constraining for 
models that explain the Fermi electron spectrum [88, 89, 9(], 72, 91]. However, this involves 



an extrapolation into the galactic center, which does involve uncertainties. As a simple 
example, we include an effect that has been heretofore ignored, namely, the variation in the 
velocity dispersion Vdisp as one approaches the galactic center. As noted earlier, many recent 



simulations with baryons have seen an increase in Vdisp roughly as r 1//4 [B3, 3C, 37] 



as one moves to the inner galaxy. As the annihilation enhancement from the long-range force 
scales as v _1 or v~ 2 in resonance regions, this can suppress the signal in the GC and GR. 

We show in figure Fig. || a comparison between a MiXDM scenario normalized to fit the 
local cosmic ray data, under the assumption that the cross section is unchanged in the inner 
galaxy, or under the assumption that it has changed from the higher velocities. As is clear 
from the figure, these effects easily allow consistency with the HESS constraints. At the same 
time, the signal in the Haze region is only affected by a factor of 1.6 or 2.4 for v^ 1 or v~ 2 , 
respectively. Of course, these curves assume extrapolation of a functional form for halo and 
velocity dispersion into regions that are not probed by simulations, and thus should be taken 
with a grain of salt 3 . 



3 In this analysis, we are scaling the cross section as (av) oc r~ 1//4 or r _1//2 . Once ira/v ~ 1, this effect 
saturates. Even in the resonance case, this does not become relevant until the inner ~ pc, and does not affect 
the plots shown. 
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Figure 7: Synchrotron radiation at 22.5 GHz from e ± of DM origin fitted to the microwave Haze 
of the same frequency. (Solid line): flux from the two XDM e species with masses 120 GeV and 
1.2 TeV, and relative annihilation rates of 1:1. (Dashed line): fit of the 1.2 TeV only. (Dotted line): 
equivalent fit for 120 GeV only. 

5. Counting the number of WIMP species with Synchrotron radiation 

One of the elements of the case for WIMP annihilation comes from the microwave excess 
in the galactic center, i.e. the "microwave Haze" |9^, ||]. Such data have already been 
used to constrain the possible models that might explain it [93]. Should Planck extend these 
data clearly to higher frequencies, it may serve to allow an additional handle on the details 
of WIMP annihilation. 

To study this, we shall primarily use the simplified model shown in Fig. |2] (XDM with 
M X1 = 120GeV and M X2 = 1.2TeV, with relative annihilations of 1:1. 

Before proceeding to the complete analysis, let us begin by studying the microwave Haze 
data we already have. As we show in Fig. ^ just using the 22.5 GHz WMAP band where 
the microwave Haze error bars are the smallest 4 , it is not possible to distinguish between the 
cases where there is one mass of M x =1.2 TeV from the MiXDM case of two 120 GeV / 1.2 
TeV. This is not surprising, as that emit synchrotron radiation at 22.5 GHz have energies 
between 5 ^ to 50 GeV, so originating from both species contribute essentially equally. 



5.1 Modeling Synchrotron Emission with Power Laws 

As we discuss in the appendix, in considering the emissivity j(v) in the range of ~ 100 GHz, 
one should expect the synchrotron emissivity Ji(v) of the full distribution of that originate 
from xi (with M Xl = 120 GeV) to be described by a steeper power law than that of \2 (with 
M X2 = 1.2 TeV), 



Ji(v) 



J 2 (v) 



~° 2 With OL\ > Ct2, 



(5.1) 



4 not including the CMB bias 
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simply because there is more energy in high energy electrons from X2 to contribute to the 
higher frequencies. That is also shown in Fig. ^ where we present the synchrotron radiation 
from the of DM origin vs latitude for v: 70, 143, 353 and 857 GHz (frequencies to be 
measured by Planck). The calculated values of the synchrotron radiation are averaged over 
longitudes of —10° to +10° as the WMAP microwave Haze is calculated from that region of 



the sky [94]. The annihilation cross sections for those particles are fixed by fitting to the 22.5 
GHz WMAP microwave Haze shown in Fig. [/]. As explained, the synchrotron emissivity of 
from a 120 GeV DM particle drops faster than that of a 1.2 TeV DM particle. 
Note that in Fig. |7| we fit the 22.5 GHz microwave Haze band data with a combination 
of the synchrotron radiation from the of DM origin and an extra constant offset, as the 
overall zero of the microwave Haze is somewhat uncertain. In Fig. |8| instead we present only 
the synchrotron radiation from the e^ 1 of DM origin, but it should be expected that with 
actual data some constant offset will be required. 



As a result of eq. 5.1 the combination distribution in the 2-species case will be described 



by a power law J comb ~ v acomb of intermediate steepness. 

Jcomb(v) ~ v- Q - mt with ai > a comb > a 2 . (5.2) 



The key point in (5.2), is that a,\ ^ a 2 . While this may appear obvious, there are actually 
subtleties involved that we shall explain. The obvious element is the direct relationship 
between synchrotron spectrum and the hardness of the electron spectrum. From synchrotron 
theory we know that for a power law distribution of electron energies of the form N(E) = 



kE p , the synchrotron emissivity of the distribution goes as [35 



J(v) oc KV (f- 1 )/ 2 . (5.3) 

Considering that for both masses we used the same annihilation channel, one would think 
that Ji(v) and J 2 (v) should be described by the same power law, which is true - up to a point. 
The injection spectra are given by the same formula^] which gives an injection power law of 0, 
and the dominant energy losses (synchrotron and ICS) scale as ^ oc E~ 2 . However, eq. ( |5.3| ) 
neglects the energy cut-off E c t in the distribution of the electrons. As we explained earlier in 



this section, it is exactly the presence of the cut-off that for frequencies v v Cc/ (5^5) where: 



3 eBE 2 c 



cf 2 2-Km'i 



sin a , (5-4) 



shouldn't be applied blindly. Thus ct\ ^ a 2 . 

A second element contributing to ot\ > a 2 is that the Thompson approximation of the 
cross-section of ICS from optical starlight breaks down at electron energies E ~ 100 GeV, 
with the actual cross section between the electron and the photon being described by the 
Klein-Nishina cross section [35|: 



3 

ctk-n = 7: a TX 



2(x + 1) 



x 2 
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Figure 8: Synchrotron radiation from the of DM origin as in Fig. ^, for 70 GHz (top left), 143 
GHz (top right), 353 GHz (bottom left) and 857 GHz (bottom right). Boost factors are from fitting 
to the 22.5 GHz microwave Haze. 



where x 



tj. At the limit where i> 1 equation 5.5 simplifies to: 



<?K-N 



\tJTX~ 1 (hx{2x) + . 



(5.6) 



This results in an energy loss rate ^ due to ICS from optical starlight that scales as E@ 
with (5 < 2 and gradually decreasing with increasing electron energy E. As a consequence, 
ICS is a less efficient energy loss for the higher energy electrons, resulting in more energy 
(fractionally) being lost to synchrotron radiation, and a harder synchrotron spectrum. 

In Fig. ^ we show the ratio of synchrotron radiation flux at frequency v averaged between 
6° to 15° in latitude and —10° to +10° in longitude (microwave Haze region), to the same 
averaged flux at 22.5 GHz. 

/ J(v)dndr 



f J(v=22.5 GRz)dndr 



(5.7) 
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Figure 9: Plot (a): Ratio of synchrotron radiation at frequency v over radiation at 22.5 GHz from 
e ± of DM origin. Assumed magnetic field model: "B-field I" (Solid line): ratio for the two XDM 
species with masses 120 GeV and 1.2 TeV, at (xi'-X?) °f 1 ; 1. (Dashed line): ratio for the 1.2 TeV only. 
(Dotted line): equivalent ratio for 120 GeV only. Plot (b): Same as for (a) but with magnetic field 
model "B-field II". Plot (c): Same as (a) but with magnetic field model "B-field III". 



where dr refers to the integration over the line of light and d£l to the integration over the 
sky region. That ratio is calculated for the cases of XDM channel with M x = 120 GeV, 
M x = 1.2 TeV separately as well as the MiXDM case. We use the 5 WMAP frequencies: 22.5, 
32.7, 40.6, 60.7 and 93.1 GHz @ and the 9 Planck frequencies 30, 44, 70 GHz (measured from 
the Low Frequency Instrument) |)7], [9^] and 100, 143, 217, 353, 545 and 857 GHz (measured 
from the High Frequency Instrument) J99|, |100| . We also show the indices of the power laws 
that fit our calculations from the lowest frequency (22.5 GHz) up to ~ 100 GHz, where 
deviations from the one power law description arise. 

Based on ( [5.2; ) if we are able to measure the power law index a of J(v), and if we know the 
mass M x (or M X2 for the 2-mass case) and the annihilation channel via which e are produced 
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then we can discriminate between the 2-species xii X2 case and the case where there is only 
X2 in nature. As shown in Fig. [9|a 5 , the power law for frequencies up to ~ 100 GeV, is for 
the case of a single species M x = 1.2 TeV equal to a = 0.40 while for the 2-species case with 
M X1 = 120 GeV and M X2 = 1.2 TeV it is a = 0.46, because of the contribution from xi that 
falls with a steeper power law (a = 0.52). Also in Fig. [9| notice that at v greater than ~100 
GHz, the power-law approximation of the distribution stops being valid. 

In general, for the MiXDM scenario in which masses are separated by an order of mag- 
nitude or greater, the difference in the power law index compared to the single WIMP cases 
is ~ 0.1 or greater. There are many uncertainties and challenges in extracting the Haze from 
future Planck data, so it remains to be seen whether this will be adequate. In this paper 
we do not intend to discuss the uncertainties arising from understanding the various astro- 
physical contributions (free-free emission, synchrotron radiation from primary electrons and 
thermal and spinning dust) ]i~0l| , 94]. Currently only the 22.5 and 32.7 GHz from WMAP 



have small enough error bars to be used in order to calculate the power law a. In |93|], the 
authors find that the ratio of the averaged Haze flux at 22.5 GHz to that at 32.7 GHz is 
J(22.5GHz)/J(32.7GHz) 1.18 ± 0.10 that gives a range, 0.2 < a < 0.6. Since pi used 



only the statistical errors of the microwave Haze, while there are additional large uncertainties 
arising predominantly from the microwave Haze's correlation with the CMB. Thus the actual 
range of values of a from the microwave Haze is even greater than this. At the moment, one- 
and two- species scenarios are indistinguishable within errors. 

Another difficulty in determining whether we are in a MiXDM scenario is the uncertainty 
in the annihilation channel, itself. Constraints from fitting the Fermi/HESS electron data 
require 1 TeV < M X2 < 4 TeV, depending on the annihilation mode, with some significantly 
softer than others. Such changes will clearly impact the spectrum of the synchrotron emission 
J(v). To explore this, we show in Fig. [K^ the ratio of synchrotron radiation flux at frequency 
v (normalized to 22.5 GHz) averaged over the microwave Haze region (as in Fig. ^) for four 
different annihilation channels for single species of mass M x . Those four channels are XDM 
e^, XDM n^ 1 , XDM 2-step cascade to and XDM 2-step cascade to fi^. The values for M x 
coincide with the values of M X2 used in Fig. ^, so give good fits to the Fermi data. 

As one can see from Fig. |l(Ji (B-field I) the power laws describing J(v) until ~ 100 GHz, 
vary between: 

0.4 < a < 0.46 . (5.8) 



The channel XDM is not shown in Fig. |§ but gives a = 0.43 (for M x = 3.1 TeV). 

Thus any channel that can be thought (with respect to the injection spectrum of the produced 
e^ 1 ) as a combination of XDM e , piF and ir^ will have a power law in the region of (|5[ 
provided x can fit the Fermi data and not violate the HESS upper limits 6 . 



See section 5.2 for definitions of magnetic fields 
6 t channels can fit the local e spectra, but are constrained from diffuse 7-rays, and we do not consider 
those modes. 
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Figure 10: Plot (a): Ratio of synchrotron radiation as in ^a (used magnetic field model: "B-field I"). 
(Solid line): ratio for the XDM e species with M x = 1.2 TeV. (Long dashed line): ratio for the XDM 
species with M x — 2.7 TeV. (Short dashed line): ratio for the XDM 2-step cascade to e ± species 
with M x = 1.6 TeV. (Dotted line): ratio for the XDM 2-step cascade to fi^ species with M x = 3.6 
TeV. Plot (b): Same as in (a) with magnetic field model "B-field II". Plot (c): Same as in (a) with 
magnetic field model "B-field III" . 



The region of ( |5.8[) is quite narrow 7 , making it a challenge to determine the mass and 
the annihilation channel based on the microwave Haze data only. Yet this information can be 
used in understanding the nature of DM. To do this, we will require an accurate measurement 
of the microwave Haze at frequencies ~ 100 GHz, which is challenging due to the various 
uncertainties involved in subtracting the astrophysical backgrounds and the CMB from the 
WMAP and Planck data. 



7 That is expected from the fact that a ~ (p — l)/2 with p the index to the spectrum of the electrons 
responsible for the radiation (see (|5.3|)) 
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5.2 Dependence of the B-field vs ISRF assumptions 



Prom the Eq. |B,1| we know that the synchrotron emissivity j'(v) of an electron inside a magnetic 
field is proportional to the amplitude of the B-field. If the microwave Haze can be explained 
from DM, then the region of the Galaxy that provides the major contribution to it is 1 < r < 2 
kpc. For that region there are significant uncertainties in the literature about the value of the 



B-fields amplitude and shape(see [102, 103, 104, |M ft0|, g07|, [l0§ and [10g] among others). 



To study the effects of the uncertainties of the B-field amplitude, we follow |74j and use 
variants of the default magnetic field profile from GALPROP: 

B(p, z) = B exp ( Rq ~ p \ exp (--) . (5.9) 



with Bq = 5pG being the value for the total(combined random and large scale "ordered") 
local value of the magnetic field, R® = 8.5kpc, p c and z c being the characteristic scales along 
p and z directions respectively. We used three options for the values of p c and z c . 

1. p c = 4.5kpc and z c = 2.0kpc, defined as "B-field I", 

2. p c = lO.Okpc and z c = 2.0kpc, defined as "B-field II" and 

3. p c = lO.Okpc and z c = 4.0kpc, defined as "B-field III". 

The magnetic field model defined as "B-field I" , provides agreement with the 408 MHz Haslam 



c 



map [110, 111], while the model defined as "B-field II", is the "conventional" galprop assump- 
tion. "B-field III" was chosen to take into consideration uncertainties in the scale height z, 
of the magnetic field 8 . In generating the plots shown in Fig.||, we used the "B-field I" model. 

In Fig. |]b and c and in Fig.lCb and c, we show the dependence of the magnetic field on 
the ratio of synchrotron radiation for the cases studied in Fig. ||a and [Toja. We can see that for 
the cases where p c = lOkpc, the power-law spread is greater than that of p c = 4.5kpc. In the 
case of p c = lOkpc, the magnetic field at 1 < z < 2 (kpc) is weaker by a factor of 2, compared 
to that of the p c = 4.5kpc case, thus the power in synchrotron radiation is by a factor of 4 
lower. As explained previously, one reason for the difference in the observed power-laws is 
the fact that ICS losses compete with synchrotron losses. Lower synchrotron radiation losses 
result in ICS losses being more important 9 . Thus the difference in the observed power-laws 
shown in and c( |lO| b and c) vs ^a ( |l^), is enhanced. 



Moreover, the weaker magnetic fields shown in |9|b and c( lCb and c) versus that shown 
in ^a ( [lO|a) results in a lower power-law values for a. This simply arises because of the 
lower total synchrotron energy losses (as we keep our ISRF constant) resulting in a smaller 
depletion of the higher energy electrons. 



8 Values of z c ~ lkpc have been used by |l02| and [ 10S ] . Such a low scale hight for the magnetic field would 
result in a greater spread of power-laws as those presented in Fig. ^ and [| and thus would make it easier to 
separate among different scenarios. 

dE 

9 dE ICS is greater 

dt total 
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As we see the current uncertainties in the B-field are significant enough to make an anal- 
ysis based solely on the observed power-laws a degenerate among the various DM models. 
The ISRF and as well the diffusion properties of the electrons are also not well modeled yet, 
adding additional uncertainties. Yet, as diffusion models improve and are better constrained, 
and as more accurate measurements of the ISRF and the B-field in the regions of interest ap- 
pear, these uncertainties will diminish. Furthermore, one would hope that in these scenarios, 
we may have some estimate of the(some) DM species from direct searches, or other indirect 
data. Thus, in concert with such results, such an analysis as we have described would be a 
complementary tool for understanding the Dark sector. 



6. 3 masses 



Of course, once two masses are present, we can consider more, and even heavier DM species. 
We will confine our discussion to the case where there are three DM species and refer to 
species %i as the lightest, \2 as the intermediate mass species and X3 as the heaviest species. 

M X1 < M X2 < M X3 . (6.1) 

As explained in section 2.1 the three species are expected to have similar fluxes of at 



E<M X1 . 

A scenario where the three particles are spaced within 100 GeV-3 TeV is essentially indis- 
tinguishable from the two-mass case. A scenario where M Xl ~ 10 GeV would be difficult to be 
seen from the flux, as Solar modulation is significant at E ^ 10 GeV, so fluxes at those 



energies can't be measured accurately. For instance, 21] considered a case of two masses 
where the lightest at 3 — 10 GeV is suggested in order to explain via its elastic scattering to 
Nal target the annual modulation signal measured DAMA/LIBRA |112| , but such a scenario 
would be difficult to extract from the signals we have discussed. We do note, however, that 
such a light WIMP could, in principle, produce an additional component of electrons that 
could contribute to the low energy Fermi Haze, although additional astrophysical sources may 
contribute as well. 

The most interesting from the point of view of locally measured fluxes is the case 
where: 

M X1 ~ lOOGeV , M X2 ~ ITeV andM X3 ~ lOTeV . (6.2) 
where additional contributions at the highest energy can arise. 



In Fig. 11 we present the case of a M Xl = 200 GeV annihilating through XDM channel, 



M X2 = 1.0 TeV XDM e^- channel as well and a M X3 =8.0 TeV annihilating through 2-step 
cascade to muons. We note that even though there is good agreement to the Fermi data, the 
projected positron fraction is below the one measured by PAMELA by almost a factor of two. 

While the two-mass case can be easily realized, the three-mass case described by ( |6.2| ) 
is actually allowed only with very specific annihilation modes. The reason behind that is 
the existence of the HESS e + + e~ upper limits H, i, that suggest a turn-off of the source 
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Figure 11: Fitted e + + e spectra as in Fig. | for annihilation of 3 XDM species with M Xl = 250 
GeV, M X2 = 1.0 TeV through e ± channel and M X3 = 8.0 TeV through 2-step cascade to /i ± . 

of responsible for the excess compared to the conventional backgrounds at ~ 1 TeV. (A 
change in the background power law at ~ ITeV would also be a reasonable possibility.) As 
a M X3 ~ 10 TeV would in general create significant fluxes of e^- at ~ 1 TeV, unless the 
annihilation mode is adequately soft, there are conflicts with the HESS upper bounds. 
Going to higher than M X3 ~ 10 TeV these tensions are aggravated, as the propagated 



contribute to the e flux with E < 5 TeV. The local primary electron flux at energies E>1 



TeV is mainly influenced by local and recent SN explosions 1 113 , 114 , 115 , 116 , 117 , 118 ], and 
very young pulsars |119| , |120| , |121] , 122 , 68, ^] as well. These uncertainties make a detection 



of a heavier mass difficult in the data. 



7. Scattering and Effects on the DM Halo 

The presence of "long-range" forces in the dark sector can be constrained from the stability of 

However in the context 



126 | 



DM halos [123, 124] and other properties of galaxies [124, \\2l 
of MiXDM there is an intriguing twist. It is possible that the heavy WIMP could have a lower 
scattering rate than the lighter WIMP. This may be because of the different mass splittings 



5, or the different mw 



im 



or, possibly, because the WIMPs couple to different force carriers. 



There are other possibilities, for instance that there are other nearby states, or that there is 
a residual symmetry allowing degenerate excited states with a long range force between them 
(analogous to the charged wino, for instance), in which case these light WIMPs might scatter 
very frequently. Such scatterings can easily be strong within the dark sector without leading 
to a significant WIMP-nucleus interaction, as was explored in | |32| |. 

In such a case, if the light WIMP scattering rate is large enough while leaving the heav- 
ier WIMP scattering rate small enough, the light WIMP halo would collapse, leading to a 
completely non-Maxwellian velocity distribution. Since this WIMP is only 1% of the total 
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energy density, however, this effect would be of little consequence for the total halo. On the 
other hand, the consequences for direct detection experiments would be profound. The num- 
ber density and velocity distribution of the lighter WIMPs locally would change significantly, 
with possibly only the highest velocity particles which are just now falling into the Milky 
Way present in the event that the light halo has collapsed. In such a case, the lightest WIMP 
might dominate the annihilation signal from the inner galaxy, while the heavier WIMP would 
dominate the local cosmic ray signals. 

It is not clear how to approach this quantitatively, except to recognize that in MiXDM 
scenarios the velocity and density structure of the subdominant species can be radically 
altered, with the subsequent consequences for direct and indirect detection experiments. 

8. Conclusions 

Cosmic ray signals from PAMELA and Fermi among other have prompted a reexamination of 
our assumptions of the nature of dark matter. Two of the most central are the assumptions 
that annihilations are into standard model states, and that there is a single parity, namely 
R-parity, that stabilizes the WIMP. 

Relaxing the latter opens the possibility that the dark matter is composed of multiple 
particles, with different masses and relic abundances. Such as scenario has been previously 
motivated to explain both DAM A and INTEGRAL. We have argued that such a "MiXDM" 
scenario, in which DM annihilates as in standard XDM into light force carriers can naturally 
lead to similar annihilation rates, even if the different WIMP relic densities are different by an 
order of magnitude. The consequence is to modify possible particle physics interpretations of 
the spectra. Softer spectra can arise naturally from the presence of the second, lighter WIMP 
state, while even bumps and dramatic spectral features can arise in the positron and electron 
spectra from the presence of these lighter states. 

We find that these models are in good agreement with the Fermi e + + e~ flux and the 
PAMELA positron fraction, with the heaviest mass restricted by the HESS upper bounds. 
These cross sections are very similar to what is needed to explain microwave and gamma-ray 
excesses in the inner galaxy region, while still below bounds from the CMB. Limits from the 
Galactic Ridge, while possibly constraining, are weakened when velocity-dependent effects 
are taken into account. 

The presence of the lighter state is possibly detectable in a variety of ways. Beyond 
the possible features in the e + and spectra, the changes in the synchrotron emission 
spectrum can also help distinguish it. In particular, synchrotron radiation from of DM 
origin at frequencies up to ~ 500 GHz, probed by Planck, are potentially different for these 
models. Unfortunately, the uncertainties in the magnetic field, ISRF and electron diffusion, 
as well as the great number of possible particle physics models, currently it is not possible 
to discriminate from synchrotron emission alone between a one DM and a two DM species 
scenarios that both fit the PAMELA and Fermi data. Yet, when combined with a better 
understanding of all the factors involved in the calculations of the microwave Haze and the 



-21- 



synchrotron radiation(B-field, ISRF, diffusion models, dust emission, Planck instrumental 
errors), such a tool may offer an important handle, especially in combination with other 
indirect and direct DM searches. 

Interesting additional possibilities exist as well. One can consider a third, heavier WIMP, 
although we find that its effects are highly constrained by the HESS electron limits, and that 
it is difficult to generate any features at high energies consistent with data. A possibility 
exists for subdominant components to have very large scattering rates, potentially collapsing 
the lighter WIMP halo, and radically altering the predictions for direct detection scattering 
rates. 
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A. DM Halo and CR propagation 



For all our calculations we considered an Einasto profile using the parametrization of [127]: 



p(r, z) 



p exp 



2 
a 



(A.l) 



with a = 0.17, r_2 = 25kpc and Rq = 8.5 kpc. We considered the local DM mass density 
to be po = 0AGeVcm~ 3 that is in agreement with calculations of the local dynamical mass 
based on Hipparcos data and the estimates on the local baryonic mass |H, |12|, |3|, |l3ll and 



also with |132j . 

While we assume that our boost from the thermal cross section is purely from particle 
physics, substructure in the distribution of DM can also exist, affecting the amplitude and 



spectrum of the signal. For instance, [133, 134, 73] find that substructure can create spectral 
hardening of the propagated spectrum at high energies. If such substructure gives a 
significant contribution to the local flux, it would also enhance the bumps at E ~ M Xl , 
although this depends on particle physics and astrophysics parameters. 

To determine the cosmic ray properties, we have calculated the propagation of the 
in the Galaxy using GALPROP|135|. In our models for the background fluxes we have 
considered a magnetic field of a local value of 5pG with a peak value of 33pG in the galactic 



center, and the standard GALPROP Interstellar Radiation Field (ISRF) [136] . For the diffu- 
sion of the CRs galprop considers a diffusion zone within which the spatial diffusion coefficient 
is considered homogeneous in space depending on rigidity as: 



Oxx 



R 

D rigid br 



D a 



(A.2) 
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with /3 = v/c, R = 2 the rigidity, Dq xx diffusion coefficient at reference rigidity D r i g idbr and 
D g diffusion index. Our assumptions for Dq xx and 5 provide good agreement with proton, 
He, C, Fe, B/C and subFe/Fe fluxes, as well as antiproton over proton ratio [|73| . For a more 
extensive discussion of background electrons and positrons one can refer to [135, 0, 137]. 



We also define the boost factor that we are going to refer to as just the overall en- 
hancement in the ^i 71 ^ "* I v I) needed to fit the data at hand compared to its value 
assuming an Einasto profile with po = 0AGeVcm~ 3 and thermally averaged cross sections 
such that the total relic mass density agrees with the current constraints from WMAP, on 



nh 2 = 0.1099 ± 0.0062 Pal 



For the case of M Xl <C M X2 and for an Einasto profile as we use in this paper, we can 
simply write the B.F. as just 10 : 



(<J2 | v |) 



B.F. « — z 1 "l /Flt t ed , , (A.3) 

M 2 

as xi contributes only~ j-fiy- to the relic mass density. In our calculations we took into 

X2 

account corrections in the relic density that scale with the logarithm of the masses M x 
B. Review of Synchrotron Emission 



The emissivity of a single electron to synchrotron radiation is given by [95]: 



v / 3e 3 Bsina v . 
J(V) = S^ocme ~ / K5/3(z)dZ (B - 1} 



v/v c 



where B is the magnetic field and a the angle between the B-field and the velocity of the 
electron, K 5 / 3 is the modified Bessel function of the second kind of order 5/3 and v c is the 
critical frequency: 

Vc= 2 7 2^: sma - (B - 2) 

which scales as cx 7 2 , as the gyroradius of the electrons in the ISM is a 7. The maximum 
value of j(v) is at v pea k = 0.286v c , at v = 3.38v c the j( v ) has 1/10 of its peak value, while at 
the very low and very high frequency limits it is given by: 



sin a 

■H v ) = ~^~2 x S 



b ^) V3 forv«v c 



v^rxi/3) V 2Vc 



(B.3) 



Considering the exponential cut-off at frequencies much higher than the critical frequency 
we can see that at the highest frequencies to be probed by Planck (up to 857 GHz), the less 
energetic electrons are not going to emit significant synchrotron radiation. For instance, for 
a ~50 GeV electron, and a magnetic field of ~ 20/iG such as expected in the Haze region, 
the critical frequency is, v c ~1 THz, thus its emissivity peaks at 300 GHz. A 1 TeV electron 
though has for the same B-field a critical frequency of v c ~400 THz thus the emissivity peaks 

at ~120 THz. 

10 within logarithmic corrections on M X2 
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